Pygmy resonances probed with electron scattering 
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Pygmy resonances in light nuclei excited in electron scattering are discussed. These 
collective modes will be explored in future electron- ion colliders such as ELISe/FAIR 
(spokesperson: Haik Simon - GSI). Response functions for direct breakup are explored 
with few-body and hydrodynamical models, including the dependence upon final state 
^ ' interactions. 

t> 

^ ' 1. Introduction 
O 

Reactions with radioactive beams have attracted great experimental and theoretical 
interest during the last two decades [ [T]. Progresses of this scientific endeavor were 
reported on measurements of nuclear sizes [ [2] , use of secondary radioactive beams to 
elucidate reactions of astrophysical interest [ El S], fusion reactions with neutron-rich 
nuclei tests of fundamental interactions [[7j, dependence of the equation of state 

of nuclear matter upon the asymmetry energy [[8], and many other research topics. 
^ ■ New research areas with nuclei far from the stability line will become possible with 
newly proposed experimental facilities. One of the projects at the future FAIR facility of 
the GSI laboratory /Germany is the study of electron scattering off unstable nuclei in an 
electron-ion collider mode [Ej. A similar proposal exists for the RIKEN laboratory facility 
in Japan [[10]. By means of elastic electron scattering, these facilities will become the main 
probe of charge distribution in unstable nuclei [[HI [12]. Coulomb excitation has been very 
useful in determining the electromagnetic response in light nuclei [d3]- For neutron-rich 
isotopes [H3j the resulting photo- neutron cross sections are characterized by a pronounced 
concentration of low-lying El strength. But it is well known that non-perturbative effects 
leading to distortion of the energy spectrum of the fragments interacting with the target 
(see, e.g. ref. [US]) is a problem of difficult nature. The nuclear response probed with 
electrons is free from such effects. 

The interpretation of the low-lying El strength in neutron-rich nuclei engendered a 
debate: are these "soft dipole modes" just a manifestation of the loosely-bound character 
of light neutron-rich nuclei, or are they the result of the excitation of a resonance? 
ITBI El [HI tlH] ■ The electromagnetic response in light nuclei, leading to their dissociation, 
is related to the nuclear physics needed in several astrophysical sites [ [31 HI [15] . The 
existence of pygmy resonances have important implications on theoretical predictions of 
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radiative neutron capture rates in the r-process nucleosynthesis, and consequently on the 
calculated elemental abundance distribution in the universe 

2. Inelastic scattering in electron-ion colliders 

Here, Jj (Jf) is the initial (final) angular momentum of the nucleus, (E, p) and 
(E', p') are the initial and final energy and momentum of the electron, and (go, q) = 
((E — E')/hc, (p — p') /h) is the energy and momentum transfer in the reaction. For low 
energy excitations, E, E' ^> hcqo, for electron energies E ~ 500 MeV and excitation 
energies AE — hcq ~ 1 — 10 MeV. 

In the plane wave Born approximation (PWBA) and using the Siegert's theorem, one 
can show that [|2T1 



da _^dN( E V(E,E 1 ,v J {EL) 



E £ IFtd " L < L) (E 7 ), (1) 



d^ldEj ^ dVtdE^ 



where a^ EL>) (E 7 ) oc dB (EL) /dE y , is the photo-nuclear cross section for the £'L-multipolarity, 
and E y = hcq . The response function, dB (EL) /dE 7 , is given by 



dB (EL) _ KJ/||yL(f)||Ji>| 



dEy 2Ji + 1 



- I drr 2+L 5$ L) (r)\ p(E 7 ), (2) 



where p (E 7 ) is the density of final states (for nuclear excitations into the continuum) 
with energy E 7 — Ef — E^. The geometric coefficient (Jf \\Y L (r)|| Jj) and the transition 
density Sp\f L ^ (r) will depend upon the nuclear model adopted. 

One can also define a differential cross section integrated over angles so that 

iJtllpEA = 2T f- M sill , dN ^f- e \ (3) 

dEy JEj/E dildEy 

and 8 m is the maximum electron scattering angle, which depends upon the experimental 
setup. Notice that the lowest limit in the above integral is 9 min = E 7 /E, and not zero. 

Eqs. UM\ show that under the conditions of the proposed electron- ion colliders, electron 
scattering will offer the same information obtained with real photons. The reaction dy- 
namics information is in the virtual photon spectrum, N^ EL ^ (E, E y , 9), while the nuclear 
response dynamics information is in eq. [2j dN^ EL ^ / dVtdE^ is interpreted as the number of 
equivalent (real) photons incident on the nucleus per unit scattering angle Q and per unit 
photon energy E~. For larger scattering angles the Coulomb multipole matrix elements 
are in general larger than the electric (EL) multipoles, and monopole transitions become 
relevant [[22]. Eq. [T]will not be valid under these conditions. 

It is found that the spectrum dN^ EIjS) (E,E 7 )/dE 7 increases rapidly with decreasing 
energies [ [21] . For E = 500 MeV and excitation energies AE = 1 MeV, the spectrum 
yields the ratios dN^/dN^ ~ 500 and dN^/dN^ ~ 100. However, although 
increases with the multipolarity L, the nuclear response decreases rapidly 
with L, and El excitations tend to dominate the cross sections. For larger electron 
energies the ratios N^ E2 ' JN^ E1 ^ and N^ m > /N^ E1 ' decrease rapidly. A comparison between 
the El virtual photon spectrum, dN^ Ei y dE y , of 1 GeV electrons with the spectrum 
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generated by 1 GeV/nucleon heavy ion projectiles was published in ref. [ [21]. The 
virtual spectrum for the electron contains more hard photons, i.e. the spectrum decreases 
slower with photon energy than the heavy ion photon spectrum. This is because, in both 
situations, the rate at which the spectrum decreases depends on the ratio of the projectile 
kinetic energy to its rest mass, E/mc 2 , which is much larger for the electron (m = m e ) 
than for the heavy ion (m = nuclear mass). 

3. Dissociation of weakly-bound systems 

3.1. One-neutron halo 

In a two-body model, the single-particle picture has been used previously to study 
Coulomb excitation of halo nuclei with success [[231 I2U [25], [26], [27] [28]. The El transition 
integral 2y . = / °° dr r 3 dpif (r) becomes 



P 2/i (S n + E r f 



., ( _M 3/2 V^(S n + 3E r 



2h 2 J -I /a x +(ir 1 E r /h 2 _ 



(4) 



in terms of the effective range expansion of the phase shift, k 2l+1 cot 5 ~ — 1/aj + r//c 2 /2. 

As shown in previous works [ [23], [21], a peak is manifest in the response function, 
dB(EL)/dE oc \l s ^ p \ 2 oc E^ +l ' 2 /(S n + E r ) 2L+2 . It appears centered at the energy [[23] 

( TP J \ 

£<P ; ~ (L + l/2)S n /(L + 3/2) for a generic electric response of multipolarity L. For El 
excitations, the peak occurs at Eq ~ 3S n /5. The second term inside brackets in eq. H]is 
a modification due to final state interactions. It is important in many cases 



3.2. Two-neutron halo 

Many weakly-bound nuclei, like 6 He or 11 Li, require a three-body treatment. In one 
particular model, the bound-state wavefunction in the center of mass system is written 
as an expansion over hyperspherical harmonics, see e.g. [[29]. For weakly-bound systems 
having no bound subsystems the hyperradial functions behave asymptotically as [ [30] 
&a (p) — ► exp (—rip) as p — > oo, where the two-nucleon separation energy is related to rj 
by — h 2 r] 2 / (2m n). (see also [E2])- The El transition matrix element is obtained by a 
sandwich of the El operator between <3? a (p) / p 5 / 2 and scattering wavefunctions. Following 
ref. [ 31J, but using distorted scattering states, 



1 ( E1 ^ = I dxdy ~/^~ y2xUp ( y ) Uq ( x ) ' ( 5 ) 

where u v (y) = jx (py) cos 5 nc — n\ (py) sin 5 nc is the core-neutron asymptotic continuum 
wavefunction, assumed to be a p-wave, and u q (x) = j (qx) cos 5 nn — n (qx) sin 5 nn is the 
neutron-neutron asymptotic continuum wavefunction, assumed to be an s-wave. 

The El strength function is proportional to the square of eq. [5] integrated over all 
momentum variables. As pointed out in ref. [ [31] , the El three-body response function 
of n Li can still be described by an expression similar to the two-body case, but with 

different factors. Explicitly, dB (El)/dE r oc EfJ^S^ + E^j . Instead of SWu one nas 

to use an effective S^J = aS2 n , with a ~ 1.5. With this approximation, the peak of the 
strength function in the three-body case is situated at about three times higher energy 
than for the two-body case. In the three-body model, the maximum is thus predicted 
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at Eq ~ 1.8S2n, which fits the experimental peak position for the n Li El strength 
function very well [[31]. It is thus apparent that the effect of three-body configurations is 
to widen and to shift the strength function dB (El) jdE to higher energies. 

3.3. The hydrodynamical model 

As with giant dipole resonances (GDR) in stable nuclei, one believes that pygmy reso- 
nances at energies close to threshold are present in halo, or neutron-rich, nuclei. This was 
proposed by Suzuki et al. [ ,33] using the hydrodynamical model for collective vibrations. 
We will use the method of Myers et al. [ 36J, considering collective vibrations in nuclei 
as an admixture of Goldhaber- Teller (GT) and Steinwedel- Jensen (SJ) modes. For light 
nuclei Goldhaber- Teller modes dominate. 

For spherically symmetric densities, the transition density, 5p p (r) = 5p p (r) Y 10 (r), can 
be calculated assuming a combination of SJ and GT distributions [I2"T]. 



where R is the mean nuclear radius, and a* is the percent displacement of the center 
of mass of neutron and proton fluids in the GT (i = 1) and SJ (i = 2) modes. For 
light, weakly-bound nuclei, it is appropriate to assume that the neutrons inside the core 
(A c , Z c ) vibrate in phase with the protons. The neutrons and protons in the core are 
tightly bound. Calling the excess nucleons by (A e , Z e ) = (A — A c , Z — Z c ), the effective 
charge for the GT mode is Z^X = (Z c A e — A c Z e ) /A. In eq. [6], ji(kr) is the spherical 
Bessel function of first order, «2 the percent displacement of the center of mass of the 
neutron and proton fluids in the SJ mode, kR = a = 2.081, and K = 2a/ jo (a) = 9.93. 

The hydrodynamical model can be further explored to obtain the energy and excitation 
strength of the collective excitations. This can be achieved by finding the eigenvalues of 
the Hamiltonian Ti. = ^aTa + ^aVa + aJ r a, where a = (ai, a^) is now a vector containing 
GT and SJ contributions to the collective motion. T and V are the kinetic and potential 
energies 2x2 matrices [ 36J. The kinetic term can be calculated from GT and SJ velocity 
fields. The potential term is related to the stiffness parameters of the liquid-drop model 
and adjusted to a best fit to the nuclear masses. The stiffness of the system is due to 
the change in symmetry energy of the system as it goes out of the equilibrium position. 
The last term in the Hamiltonian is the Rayleigh dissipation term, which is related to the 
Fermi velocity of the nucleons [ 36J and yields the width of the eigenstate. 

As shown by Myers et al. [EE], the liquid drop model predicts an equal admixture of 
SJ+GT oscillation modes for large nuclei. The contribution of the SJ oscillation mode 
decreases with decreasing mass number, i.e. a — > (ai, 0) as A — > 0. This is even more 
probable in the case of halo nuclei, where a special type of GT mode (oscillations of the 
core against the halo nucleons) is likely to be dominant. For this special collective motion 
an approach different than those used in refs. [ EH] and [ [33] has to be considered. 

It is easy to make changes in the original Goldhaber and Teller [[31] formula to obtain 
the energy of the collective vibrations, yielding 




(6) 




(7) 
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where A r = A c (A — A c ) j A and a is the length within which the interaction between a 
neutron and a nucleus changes from a zero- value outside the nucleus to a high value inside, 
i.e. a is the width of the nuclear surface, ip is the energy needed to extract one neutron 
from the proton environment. Goldhaber and Teller [[34] argued that in a heavy stable 
nucleus cp is not the binding energy of the nucleus, but the part of the potential energy 
due neutron-proton asymmetry. In the case of weakly-bound nuclei this picture changes 
and it is more reasonable to associate (p to the separation energy of the valence neutrons, 
S. I will use ip = j3S, with a parameter (3 which is expected to be of order of one. Since 
for halo nuclei the product aR is proportional to S~ l , we obtain the proportionality 
Epr oc S. Using eq. [7] for n Li , with a = 1 fm, R = 3 fm and <p = S2n = 0.3 MeV, 
we get Epr = 1.3 MeV. Considering that the pygmy resonance will most probably decay 
by particle emission, one gets E r ~ 1 MeV for the kinetic energy of the fragments. This 
is about a factor 2 larger than what is obtained in a numerical calculation [ [2T] . But 
it is within the right ballpark. The hydrodynamical model is very unlikely to be an 
accurate model for light, loosely-bound, nuclei and is significant only in that a reasonable 
magnitude of the resonance energy is found. 

Both the direct dissociation model and the hydrodynamical model yield a bump in 
the response function with position proportional to S, the valence nucleon(s) separation 
energy. In the direct dissociation model the width of the response function depends on the 
separation energy and on the nature of the model, i.e. a two- or a three-body model. In 
the case of the pygmy resonance, this question is completely open. The hydrodynamical 
model predicts [[36] for the width of the collective mode T = hv/R, where v is the average 
velocity of the nucleons inside the nucleus. This relation can be derived by assuming that 
the collective vibration is damped by the incoherent collisions of the nucleons with the 
walls of the nuclear potential well during the vibration cycles. This approach mimics that 
used in the kinetic theory of gases for calculating the energy transfer of a moving piston 
to gas molecules in a container. Using v = 3yp/4, where vp = \jlEp /m^ is the Fermi 
velocity, with Ep = 35 MeV and R = 6 fm, one gets T ~ 6 MeV. This is the typical energy 
width a giant dipole resonance state in a heavy nucleus. But in the case of neutron-rich 
light nuclei v is not well defined. But the piston model is not able to describe the width 
of the response function properly. Microscopic models, such as those based on random 
phase approximation (RPA) calculations, are necessary to tackle this problem. The halo 
nucleons have to be treated in an special way to get the response at the right energy 
position, and with approximately the right width. Right now, the problem remains if the 
experimentally observed peak in dB/dE is due to a direct transition to the continuum, 
weighted by the phase space of the fragments, or if it proceeds sequentially via a soft 
dipole collective state. 
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